OCHA-PP-277 
October 2008 



Factorization of number into prime 
numbers viewed as decay of particle into 
elementary particles conserving energy 

Akio SUGAMOTO 

Department of Physics, Ochanomizu University, Tokyo 112-8610, Japan 

Abstract 

Number theory is considered, by proposing quantum mechanical models and 
string-like models at zero and finite temperatures, where the factorization of 
number into prime numbers is viewed as the decay of particle into elementary 
particles conserving energy. In these models, energy of a particle labeled by 
an integer n is assumed or derived to being proportional to In n. 
The one-loop vacuum amplitudes, the free energies and the partition functions 
at finite temperature of the string-like models are estimated and compared 
with the zeta functions. The SL{2, Z) modular symmetry, being manifest in 
the free energies is broken down to the additive symmetry of integers, Z_|_, after 
interactions are turned on. In the dynamical model existing behind the zeta 
function, prepared are the fields labeled by prime numbers. 
On the other hand the fields in our models are labeled, not by prime numbers 
but by integers. Nevertheless, we can understand whether a number is prime 
or not prime by the decay rate, namely by the corresponding particle can 
decay or can not decay through interactions conserving energy. Among the 
models proposed, the supersymmetric string-like model has the merit of that 
the zero point energies are cancelled and the energy levels may be stable against 
radiative corrections. 



1 Introduction 



Recently, the long-standing unsolved conjectures of Fermat conjecture, Sliimura- 
Taniyama conjecture and Sato- Tate conjecture are successively solved. (See |1J, 
^ and the references cited therein.) The lecture by Kato [T] and the reading of 
the original paper by Riemann [3] are enough to move the author strongly towards 
considering the number theory, even though he is completely ignorant of the field. 
Motivation of this paper is as follows. 

In number theory, how a number factorizes into prime numbers is a key issue, 
while in particle physics how a particle decays into elementary particles is also a 
key issue. These two key issues are intimately related, if we identify the energy 
En of a particle labeled by a positive integer n = 1, 2, 3, ... is proportional to Inn. 
Namely, 

En = Ljlnn, (1) 

where a; is a unit of energy. 

Then, factorization of a number into prime numbers {pi,p2,P3, . . .}, 

n = {pir{P2r{P3r---, (2) 

can be viewed as the energy conservation law, 

En = Ci-Epi + C2-Ep2 + C3-E'p3 + • • • ) (3) 

where ci, C2, C3, . . . are positive integers. 

In this sense, a particle (labeled by) n can be a bound state made up of ci 
particle pi, C2 particle p2, C3 particle ^3, and so on, and it can decay into elementary 
particles with smaller energies. In this decay process the total energy should be 
conserved. Therefore, it may be possible to understand number theory in terms 
of (particle) physics terminologies. 

Even though the factorization of number and the particle's decay following en- 
ergy conservation law are nothing but the paraphrasing of each other, the physical 
intuitions on energy, free energy, decay width, unitarity and the optical theorem, 
supersymmetry, etc. may be useful in number theory. Also the physical treatment 
on the symmetry and its breaking of the dynamical system may be helpful there. 

Purpose of this paper is to attempt for rewriting number theory in terms of 
the quantum mechanical theory and the two dimensional string-like quantum field 
theory, at zero as well as finite temperatures. 
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2 Quantum model with infinitely many species 
of particles 



In quantum theory, one species of particle, labeled by n (= 1, 2, 3, . . .), is described 
by a field 0„(t) depending on time t. If the Lagrangian L is given for a infinite set 
of fields {0n(t)} (n=l,2,3,. . . ), the dynamics between particles, or how (with what 
probablity) each particle decays into other particles is known. The simplest quan- 
tum model "Model A" is the non-relativistic bosonic theory with cubic coupling 
A, whose Lagrangian is 

La = Y.{^nit)^J^Ut)-E^4>nmn{t)^ 

ni,n2,n3 

Here, the "overhne" means the complex conjugation and the affix indicates 
"bosonic". The canonical momentum n^^(i) of the field 0n(O is 





SO that the canonical commutation relations read 

[Mt)Am{t)] = Mt)Amit)]=0, (6) 
[<f>n{t),'$m{t)]^ Km- (7) 

If the interaction is switched off (setting A = 0, or moving to the interaction 
representation), the equation of motion determines the time dependence of the 
field explicitly, 

(j)n{t) tTSt e-'^'-'an, (8) 

with 

[an, am] = [an, am] = 0, (9) 

[an,am\^Sm,n- (10) 

From these expressions, a„ and a„ are, respectively, understood as the creation 
and annihilation operators of the particle of the species n. 
The Hamiltonian of Model A is 

ni,n2,n3 
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Therefore, the free Hamiltonian becomes 

=J2e'^ Mt)Ut) = Y.K anan. (12) 

n n 

The energy of the state |ci, C2, C3, . . .), which is obtained from the vacuum |0) by 
creating ci particle rii, particle ?t,2, C3 particle ri.3, . . ., 

|ci, C2, C3, . . y = l-(a„J^^^(a„J^^^(a„3)^^ ■ ■ ■ |0), (13) 
n\\ ^3! 

is given by 

H^^^\c,,C2,cs,...)' = E'^\c^,C2,cs,...)\ (14) 

where 

E,^ = Ci<+C2<+C3< + ---. (15) 

Here, the vacuum |0) is defined by a„|0) = for any n = 1, 2, 3, . . ., and the affix 
"6" means "bosonic" as before. This bosonic state corresponds to the factorization 
of a number n into 

n = {mrin2rin,r---. (16) 

If we can restrict the labels of the bosonic fields to be prime numbers, then Eq. f fT6l) 
reproduces precisely Eq.([2]). 

Next, we introduce the fermionic field ipn{t) labeled by a natural number n 
which interacts with the bosonic field (j)n{t) by Yukawa coupling y. Then, we have 
"Model B" whose Lagrangian is given by 



+ T.[^nit)^f^Mt)-E^i^^it)Mt)) 



^ A(0„^(t)0„,(t)0„3(t)+0„,(t)0„^(t)0„3(t)). (17) 



Here, the affix "/" means "fermionic". 
The Hamiltonian of Model B is 



Hb = Y.{K.{<Pnit)Mt))+El{iJ^{t)Mt 

n 

+ E l/(^niWV^n2(t)) (^naW+^nsW 
ni,n2,n3 

+ E A(0,^(t)0„,(t)0„3(t)+0„,(t)0,^(t)0„3(t 



ni,n2,n3 
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and its free part reads 

= E {^nit)Mt)) + El (?Jt)^n(t))) . (19) 

n 

The quantization of the fermionic field, ipn{t), is done similarly as the bosonic 
field, but anti-commutator is used rather than the commutator in the former case. 
Since the conjugate momentum of ipn{t) is 

^^-{t) = ,x = #nW, (20) 

the canonical commutation relations read 

{Mt),i^mm = Hnit),i^mm = 0, (21) 
{Mt),i^mit)} = 5n,m. (22) 

If the interaction is switched off (or in the interaction representation), we have 

Vn(t) e-'^"*6„, (23) 

with 

{bn,bm} = {bn,bm} = 0, (24) 
{bn, &m} — ^m,n- (25) 

From the anti-commutation relations, creation of fermionic particle can not be 
doubled, so that the state with the fermionic excitation takes the following form 

\di, d2, 4, . . y = bnj)n^nz " " " |0). (26) 

The energy of the state (eigen- value for H^) is 

Ei = <+E4+E4 + ---. (27) 

Therefore, the fermionic state represents the natural number n decomposed into 

n = nin2n^ ■ ■ ■ . (28) 

If we restrict the labels of the ferimonic fields to be prime numbers, then Eq. fl27|) 
reproduces the number n, which is not divisible by any square of prime numbers 
other than 1. For this kind of number n, Mobius function ^{n) is introduced in 
number theory in order to count the number of prime factors in n. This function is 
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called fermion number F in physics, representing the number of fermions excited 
from the vacuum, where each fermion can not be excited repeatedly. Then, (— 1)^ 
is the fermionic parity (—1)^ in physics. Number theory seems to be coming and 
going between the bosonic representaion and the fermionic representation (see p. 8 
of [3J), so that the "supersymmetry (SUSY)", the symmetry exchanging bosons 
and fermions {4>{t) ^ "ipit)) and is the target of the next collider experiments such 
as LHC and ILC, seems to play an important role also in number theory. 

In order to have the supersymmetric (SUSY) model "Model C" , we may tune 
the parameters in Model B as follows: 

= El (or cu^ =uj^ = u;), and A = y. (29) 

In quantum mechanics or in quantum field theory, one state changes to another 
state (making transition) in the lapse of time. The conservation of energy is not 
guaranteed at finite lapse of time because of uncertainty principle in quantum 
theory, but it holds at infinite lapse of time. Therefore, we have to consider the 
transition at infinite lapse of time also in number theory to gurantee the energy 
conservation. The transition probabitlity from an initial sate \i) to an final state 
I/) at infinite lapse of time is called S-matrix element and is given by 

{f\S\t) = 6f, + z{f\T\t), (30) 



{f\m = {f\H(%)+ Y: im^'^rn) ^ (m|gW|z) 

+ E (/|g^^Vi) ^(0) ' (0) ^ . {rm\H^''\m2) Lo) ^ . {m^lH^^'l^) 

+ (31) 

where H^^^ is the interaction Hamiltonian, the remaining part of the Hamiltonian 
other than the free Hamiltonian H^^\ which depends on the coupling constants A 
and y, while E^^ is the energy eigen- value of the state \a) for the free Hamiltonian 

The energy eigen-value (energy level) is usually shifted by the interactions, 

i?„ = + (n|g(^)|n) + (n|g(^)|m) ^ , (m|gW|n) 

Er - Ein' + le 

+ (32) 
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However, in order to study how a natural number factorizes into other (prime) 
numbers, the shift of energy level from the original value E!;^^ = lu Inn is unfavor- 
able. In this respect the supersymmetric model such as Model C is promissing, 
since in the supersymmetric model the shift of energy level does not usually occur 
owing to the cancellation between bosonic and fermionic intermediate states. 

Now, we will discuss the decay of particle. Since S-matrix is unitary 

J2{f\S^\m){m\S\z) = 6f,, (33) 

m 

the so-called optical theorem holds, 

rtotai(ri) = Y.\HT\n)\^ = 2Im{n\T\n), (34) 

m 

where the "total decay width" Ttotai{n) means the probablity of the particle n 
decaying into all kinds of final sets of particles, and is itentical to the imaginary 
part of the forward {i.e. initial and final state are identical) scattering amplitude. 

Therefore, whether a natural number n is prime number or not can be deter- 
mined by the analytic property of the forward scattering amplitude as follows: 

Im{n\T\n) = 0, if n is prime, (35) 
7^ 0, if n is not prime. (36) 

Hence the number of prime numbers that are smaller than a given quantity x, 
that is F{x) in can be given as 

[x] 

F{x) = ^6lJrn{n\T\n)fl}, (37) 
n=l 

where [x] is the Gauss symbol, giving the maximum integer less than x. The 
expression may be approximated for large 

F{x) ^ J^y 6{Im{y\T\y)). (38) 

Therefore, if the proper expression of the decay width or the forward scatter- 
ing amplitudes is obtained, perturbatively or non-perturbatively, the distribution 
function of prime numbers, F{x), can be estimated. 

3 String-like model giving logarithmic series of 
energies 

In this section we will propose string-like models, giving automatically the loga- 
rithmic series of energy. En = u In n. 
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The simplest example is "Model A'" , which reproduces Model A of the last 
section. Its Lagrangian density is 

- d - f d\ 

jCa' = (p{t,a)i—(f>{t,a)-(j){t,a)u\n[l-i—\(f){t,a) 

- {27rf /^X5{a) (^{t, a)^{t, a)^{t, a) + a)^{t, a)4>{t, a)) , (39) 

where the bosonic field a) is assumed to satisfy the periodic boundary condi- 
tion, 

0(t,(7) = 0(t,(7 + 27r). (40) 

So the physical object described by this model is a string-like one of the closed 
type, being parametrized by a. 

Using the Fourier expansion of (f){t, a) 

OO -1 

</'(^'^) = E^'^n+i(i)e+'"^ (41) 

n=0 V^TT 

the Lagrangian of the Model A' reproduces the Model A, that is, 

/■27r 

/ dajCA'{t,a)^LA{t). (42) 
Jo 

Therefore, Model A' is the field theoretical realization of the Model A. 

Here, we have to comment on what kind of set of fields has to be prepared. 
In the Model A', the set of fields prepared is labeled by the whole integers, n — 

1, 2, 3, 4, 5, 6, It is of course better to prepare the set of fields labeled by the 

whole prime numbers, p = 1, 2, 3, 5, 7, . . . for the study of number theory. To 
prepare the latter set in physics, however, we have to consider the more elaborate 
dynamical system, such as Chaos e.t.c. On the other hand, the set of fields labeled 
by the natural integers is easier to prepare, but even with this set of fields prime 
or not prime can be understood as the elementary nature of particle, or as decay 
or not decay of particles, so that we will adopt this simple set of fields, labeled by 
the natural integers in this paper. 

Next, we will reproduce the Model S by a string-hke "Model S'". Its La- 
grangian density is given by 

Cb' = ^{t,a)i^^(l){t,a)-^{t,a)uHn(^l-i-^^(l){t,a) 

— d — f d \ 

+ ilj{t,a)i—ilj{t,a)-ilj{t,a)unn\l-i—\ ij{t,a) 

- {2n f^y5{a) (^(t, a)^(t, a)) (0(t, a) + 0(t, a)) 

- (27r)3/2A5(a) (0(t, a)0(t, a)0(t, a) + 0(t, a)0(t, a)0(t, a)) , (43) 
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where the fermionic field ip{t, a) is also assumed to be periodiqj, 

^it,a) =ijit, a + 2Tr). (44) 
Similarly, Fourier expansion of a) 

OO -| 

V;(t,a) = ^-=V;„+l(t)e+^'"^ (45) 

n=0 V 

reproduces the Lagrangian of the Model B as 

r daCB'{t,a) = Lsit). (46) 
Jo 

The supersymmetric (SUSY) string-like "Model C" can be also introduced by 
tuning the parameters in Model B' as in Eq. fl29l) . 

4 One-loop vacuum amplitude, free energy and 
zeta function 

Scattering amplitudes can be estimated using the old-fashioned perturbation ex- 
pressed in Eg. (13 II) . but are obtained more easily by the path integration method 
of Feynman. For example, the partition function Zji[J] and the connected gen- 
erating functional at zero temperature (T = 0) of Model A' are given by 
introducing the external source J as 

Za[J] = e*^^'-"] = / Vmi^m exp ^ J dtdcr [Ca' + 7(00(0 + , (47) 

Here, (0 denotes (t, a). Successive functional differentiation by sources determines 
the connected N-point function (amplitude) for 0(0; 



(0(^1), 0(6),..., 0(e^)) 



(4J 

j=o 



«^5J(ei)5J(6)---^^(eJv) 

This N-point function is directly connected to the N-point scattering amplitude by 
the so-called reduction formula of Lehmann, Symanzik, and Zimmermann (LSZ). 

First, we examine the one-loop connected vacuum [N = 0) amplitude at T = 0, 
since for this case we need not consider the interactions as well as the external 
sources. 

db 

In / 

{t,a)=(f>{t+b,CT) 



iWA{l-loop vac)T=o = / V^{^)V(j){^)expi dtdaCA'- (49) 

Jo b J(f>(t,a)=(f>(t+b.a) J 



"Anti-periodic boundary condition is also possible for the fermionic field. 
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Here, one-loop means the bosonic (fermionic) field propagates and comes back 
(with opposite sign) to its original value after the lapse of time b, forming a loop, 
and so we impose the condition, 

0(i,(T) =0(^ + 6,(7), (50) 

in addition to the periodic boundary condition for the closed string-like object, 

(f){t,a) ^(f){t, a + 2Tr). (51) 

Then, the parameter space {t, a) becomes a torus, having the modular symmetry 
of SL{2,Z). 

Using the mode expansion of the field 

<l>it,<^)^ E E^«-,n+ie-'^^"^*""^^> (52) 



m=— oo n=0 

and performing the Gaussian integration, we obtain 



iiy^(l-loop vac)T=o = / T" E E ~ ( "T"^ ~ a;ln(l -|- n 



(53) 



dh 

m=—oo n=Q 

where ln(det A)~^ — —TrlnA has been used. Summation over m is carried out by 
J2 ln(m + X)= dy =ln(sin7rX), (54) 

m=-oo m=-oo ^ + 11 

where the irrelevant constants are ignored, then we have 



°° foo db ( (buj 
iiyA(l-loop vac)r=o = - E / ^ In sin -— In 

„=i -^0 b V V 2 



(55) 



The free energy (Helmholtz free energy) at finite temperature is almost similar 
quantity as the one-loop vacuum amplitude at zero temperature. At finite tem- 
perature vacuum is replaced by the heat bath. Usually this replacement is done 
by making time t to be imaginary {t — it), and the periodicity is imposed on the 
interval < r < /3, in the similar manner as the one-loop case. Here, we have to 
fix /5, corresponding to fixing the temperature, and do not integrate over (5. This 
13 can be identified to the inverse temperature (more precisely, h/ (3 — ksT, 
when recovering Planck constant/27r = h and Boltzmann constant^/ce). The va- 
lidity of this replacement can be understood from that the quantum mechanical 
time evolution operator becomes the Boltzmann weight in the finite temperature 
statistical mechanics, namely. 



e 
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Then, the free energy Fa of the Model A' is given by the partition function as 

Now, we can obtain the free energies at finite temperature of our Models A'^'\ 
B^'\ and C^'), in the case without interactions, as follows: 



Fa{I3) = ^In ^sinh ^^Inn 



°° f 1 1 
y \-ujHnn + -\n(l-e~ 



/3a;'' In n 



(57) 



FsiP) = F^(/?)-if:in(^cosh(^^lnn^^. 

^ |:{_(.'_./),„„+_,„(__^)|. (58) 

^-cW = -i:iu(^^^]. (59) 
The corresponding partition functions of Models A', B', and C are 

^^(/^) = n -"^'^^T^lr^' (60) 

72 — positive integers 

zbw = n "-""'"'--''l^^. (61) 

71 — positive integers 

1 _(_ r)-/^'^ 

^-(/^) = n (62) 

71 — positive integers 

The zeta function is a kind of partition function, being expressed as 



C(-)= E n-'= n (r4^1' (63) 

72 — positive integers p — prime numbers 



1 — p 



so that the parameter s be (3uj = uo/T in physics, namely the energy unit of 
the system huj measured in the energy unit of the heat bath ksT. Therefore, 
the discrete symmetry between s ^ (1 — s), studied by Riemann [3] means the 
symmetry between the hot heat bath {ksT > hu) and the cold heat bath {ksT < 
hu). The symmetry is established by using Jacobi's identity, or in physics, by using 
the dual transformation from the model defined on an integer lattice to the model 
defined on its dual lattice, the physical meaning of this symmetry is expected to 
be clarified. The dynamical system behind the zeta function is similar to Model 
A, but there is a difference. Contrary to Model A, the zeta function has no zero 
point energy (without having the factor n~^^) and the set of fields are labeled, not 
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by integers but by prime numbers. The reason why — iln(^(s) plays an important 
role in [3] is that it is the free energy of the dynamical system behind the zeta. 
The zeta-hke function relevant to Model B is 

c^(.)= y: n-'= n i^+p-')^ (64) 

n — fcrmionic numbers p — prime numbers 

where the affix "f" means "fermionic" and the "fermionic number" means the 
number not divisible by the square of prime numbers. The corresponding zeta 
function for Model C is ({s)('^{s). 

The partition function Z{(3) and the free energy -F(/3) are in general 

Z{f]) = Tr e-^, F{[3) = In Z{(3), (65) 

where H is the Hamiltonian operator of the system. If we replace the trace by the 
supersymmetric trace Str, then we have 

Z{(3) = Str e"^ = Tr {-ife-'', F{(3) = In Z{(3) (66) 



The example of this is 

1 



as) 



E {-lYn-% (67) 



positive integers 



and for the supersjmimetric (SUSY) model C, Zc{(3) = 1 and Fc{P) = becomes 
trivial. This means the partition function and the free energy defined by Str is 
a good measure to observe SUSY breaking. This Str appears also in the one- 
loop amplitude at zero temperature, since each particle contributes to the loop 
amplitude with the fermionic parity (—1)^. Now, the one-loop vacuum amplitude 
for the Model B and C reads 

^ r°° dh /sin Inn) \ 
zWb{1-Ioop vac)T=o = - E / T ■ , ' (68) 

lo b \ sin Inn)/ 



n=l 



2 



iWcil-hop vac)T=o = 0, (69) 

where it it noted that (—1)'^ alters "cos" to "sin". Now, the vanishing of the one- 
loop vacuum amplitude in Model C, Wc = 0, can be shown, which is, however, 
an entrance to the general statement in SUSY model that there is no radiative 
corrections to the quadratic mass terms of the fields. This general statement is 
expected to hold also in our Model C with the proof. 
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5 Symmetry of the dynamical system 



One of the merit of introducing the string-hke models (Model A', Model B', Model 
C", etc.) is to elucidate the symmetry of the dynamical system, such as SL{2, Z), 
which is implicit and hidden in the quantum mechanical models (Model A, Model 
B, Model C, etc.), even though the latters are equivalent to the formers. 

As was stated in the last section, the symmetry of parameter (coordinate) space 
(r, cr) which appears in the one-loop vacuum amplitues or the free energies of the 
string-like models is SL{2, Z), the symmetry of the torus, 

where the matrix T 

is an element of SL{2, Z). Corresoponding to the coordinate transformation, we 
have 

where (6, t) are understood to be (i/3, ir) for the free energies. 

If number theory is considered as a dynamical system in physics, its symmetry 
appears in the action or in the free energy. Under SL{2, Z) transformation of 
coordinates, the fields are assumed to transform as scalars, 

0(0 ^ 0'(o = m, and m ^ V''(o = m- m 

This is the ordinary treatment of the coordinate transformation such as the general 
coordinate transformation (the differomorphism) in physics. Then, the free energy 
of the Molel A', for example, transforms as 

Fa{(3) - Fa{[3)' = ^ E 1^ (f ^' - ln(l + n')) , (74) 
where (m', —n') are given by 
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For the S'L(2, Z) transformation, if J2m,n can be replaced by J2m',n'y then the 
invariance of the free energies under SL{2, Z) hold. This is the case when the 
interactions are switched off (or A = and y = 0). 

In order for the replacement of J2m,n by J2m',n' works correctly, the range of m 
and n should be both — oo < 171,11 < +00. So far the range of is < n < +00, 
the modification of the model is necessary. This can be done naively by using 
ln(l + n) = ln(l + |n|), but the more careful examination is needed based on the 
analytic continuation of the logarithm. 

When the interactions are turned on, however, the symmetry is broken down to 
the additive group of integers Z+ which is defined by imposing 7 = and \6\ = 1. 
The reason is as follows: The interaction terms have 5(cr) which is not invariant 
under SL{2,Z), but becomes invariant, 

5ia') = 5{27c-ft/b + 5a) = ^<5((t) = 5((t). (76) 

when 7 = and \6\ = 1. This shows the dynamical system is invariant under Z_|_ 
even after the interactions are turned on. 

In order to understand SL{2, Z) in the same manner as the discrete symmetries 
such as the parity (P), the charge conjugation (C), and the time reversal (T) , the 
following treatment is convenient. This is to transform only the fields, without 
changing the coordinates. Namely, SL{2, Z) transformation is carried out by 

^{T,a) ^ 0(r',a'), (77) 
V^(r,(T) ^ ^(r',a'), (78) 

where the (t', a') and (t, cr) are related by Eq. (!70l) . This treatment is, of course, 
equivalent to the former. 

Then, the tranformation can be understood as the transformation of m and n, 

SL{2, Z) : (m, -n) (m', -n') = {in, -n)T. (79) 

Now, the discrete symmetries of P, C, and T are 

P : (0 and i)){T, a) {<p and V)(t', a') = {(p and ^)(r, -a), (80) 
C : (0 and ^)(r,a) ^ (0 and V)(r',0 = (0 and ^)(r, a), (81) 
T : (0 and ^) {r,a) ^ (0 and iP){t', a') = (0 and ^){-t,cx), (82) 
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which can be understood as 

P : (m, n) — > (m', n') = {m,—n), (83) 
C : (m, n) — > {m',n') = {—m,—n), (84) 
T : {m,n) {m',n) = {—m,n). (85) 

Then, the extension of the range of n to —oo < n < +oo makes the dynamical 
system to be invariant also under P and C. so that the concept of the "anti- 
paritcle" labeled by —n, the charge conjugate state of the particle labeled by n, 
appears here. The energy levels of particle and anti-particle are usually identical. 
This implies the above mentioned naive ln(l + n) = ln(l + |n|), but the analytic 
continuation often requires the additional iivr, making the imaginary part of the 
energy to be asymmetric, or introducing the unusual structure in the Fermi (as 
well as the Bose) surface. 

6 Conclusion 

In this paper, factorization of a number into prime numbers is considered as a 
particle decays into elementary particles conserving energy. For this purpose, 
wc have to prepare infinite spicccs of particles having energy En = ujlnn {n = 
1, 2, 3, 4, ...). These infinite spicccs of particles arc combined to a single field, giving 
the stirng-like models. By introducing the proper interactions, wc can understand 
whether a number is prime or not by the fact that the corresponding particle can 
or can not decay. This is a story at zero temperature. 

After moving to finite temperature, and estimating the partition functions and 
the (Helmholtz) free energies, we can understand the (Riemann) zeta function is 
a kind of partition function for the dynamical model consisting of a set of fields 
labeled by prime numbers. The models proposed in this paper, however, consist 
of a set of fields labeled, not by prime numbers but by integers. 

Models made of fields labeled by prime numbers are more direct in the context 
of primencss. On the other hand the models proposed in this paper arc made of 
fields labeled by integers, so that they are less direct in this primeness. However, 
the models are field theoretically constructed and the interactions are introduced 
naturally to study the transition (or decay) probablity of particles, so that the 
deficit above may be covered. For example, the distribution function of prime 
numbers can be estimated, if the decay widths or the forward scattering amplitudes 
of particles are known. 
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The zeta function is a kind of finite temperature partition function of the 
dynamical system. We know the finite temperature statistical mechanics is a 
powerful method to study the averaged or the smeared property of the system 
surrounded by the heat bath. By the existence of the heat bath, the energy 
conservation is violated, since the heat bath emits or absorbs energy. Therefore, 
in additon to this averaged or smeared study of number theory using the free 
energies or the zeta functions, we wish to utilize the zero temperature field theory 
in which the energy is always conserved and the primeness is also recognized by the 
transition probability. Therefore, going back and forth between zero and non-zero 
temperature studies is necessary. In these studies, the supersymmetric (SUSY) 
models such as Model C", may be interesting, since the supersymmetric model has 
the merit on the stability of energy levels against radiative corrections and the 
cancellation of zero point energies. 
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